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I. INTRODUCTION 



In feedback control systems, the plant control is often 
a function of the current plant state vector. A knowledge 
of the state vector is therefore required for the controller. 
In most systems, however, the entire state vector is not 
available through direct measurement. It is therefore 
necessary to estimate the unknown states. An observer is a 
device which can be used to reconstruct all or part of the 
state vector [1] . 

Traditionally, observers and controllers have been 
designed separately. Time and frequency domain specifica- 
tions are satisfied through the selection of controller 
feedback gains as if all of the plant states are available. 
The observer is designed to estimate the plant states by 
proper selection of the observer gain matrix. The observer 
output is used to supply estimates of the missing states to 
the controller [2] . 

In this thesis, after presenting the necessary background 
material, a method for designing observers which minimizes 
a quadratic performance measure is presented. The method is 
also extended to the simultaneous design of the observer and 
the controller. For simplicity, only linear, time- invariant 
systems are considered. The design method relies on the use 
of a digital computer algorithm, however, the computations 
are performed off-line. Numerical examples are presented 
in order to illustrate the new technique. 
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II. MATHEMATICAL BACKGROUND 



A. TRADITIONAL DESIGN TECHNIQUES 

th. 

The dynamic characteristics of. an n -order, linear, 
time-invariant, free system can be described by the linear, 
first order differential equations 

£(t) = Hs(t) (1) 

where s_(t) is an nxl vector and H is a constant nxn matrix. 
The response of the system to any non-zero initial condition 
is determined by the elements of the matrix H which are, in 
turn, a function of the system configuration and physical 
constants. By adjusting system parameters, a suitable H ma- 
trix is realized which causes the system to meet certain 
system performance specifications. 

Traditionally, performance is specified as a consistent 
set of time and frequency domain criteria. Typical speci- 
fications are given in terms of rise time, settling time, 
percent overshoot, steady-state error, gain margin, phase 
margin, and bandwidth. These criteria are then translated 
into a control strategy through the use of some relatively 
simple mathematical tools; generally by trial-and-error 
techniques. Appropriate assumptions are usually made to 
attempt to approximate complex, higher-order systems by 
simpler, lower-order models that can be handled with greater 
ease . 

This approach, although empirical in nature, and rely- 
ing heavily on the ingenuity and system know-how of the 
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designer, is relatively straightforward for simple, single- 
input, single -output systems. Many modern systems are 
highly complex, however, with multiple inputs and multiple 
outputs. These systems do not readily lend themselves to 
the mathematical techniques alluded to earlier and often 
performance is not easily specified in terms of the tradi- 
tional criteria. It is therefore desirable to have an al- 
ternative method for prescribing desired system performance. 

B. THE QUADRATIC PERFORMANCE MEASURE 

One such index of system performance is the quadratic 
performance measure 



which assigns a positive scalar value to the infinite time 
integral of the time -multiplied quadratic form in s^(t) and 
Q. The elements of the nxn matrix are specified to deter- 
mine the relative contributions to J of the elements of the 
vector s_(t) . The value of the time-multiplier exponent, r, 
establishes the relative contributions to J of the earlier 
and later values of the elements of s_(t) . Interpreting J 
as a cost associated with a particular system configuration, 
it becomes the task of the designer to find values for the 
elements of H which are at his disposal that yield the lowest 
value for J. This choice for H is said to be optimal since 
it minimizes the cost. 

Examples of systems to which a performance measure of 
the form of (2) can be applied arc the observer error state 
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( 2 ) 



o 
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equations 



£(t) = Fe(t) 



(3) 



and th.e state-error differential equations 
plant -observer /controller 



x(t) 




A+BK | BK 




x(t) 


e(t) 




0 | F 




e(t) 













of a combined 



(4) 



Equations (3) and (4) will be developed in detail in Sections 
III and IV, respectively. 



C. EVALUATING THE QUADRATIC PERFORMANCE MEASURE 

In Ref. 3 it is shown that if (1) is a stable,'*' linear 
system and if Q is any positive definite matrix, the per- 
formance measure (2) can be expressed as 



J = (-l) r+1 r:s T (0)V r+1 s(0). (5) 



The real symmetric matrix V_ r+ -^ is the unique solution of 
the simultaneous linear, algebraic matrix equations 



H V Al + V ,H = V 
s + 1 — s + 1— — s 



s = 1 , 2 , . . . r 



,T 



( 6 ) 



H V x + V^H = -Q . 



It can be shown [4] that if H is a stability matrix and if 
Q is positive definite then the matrix V^ is also positive 
definite. The argument can easily be extended to show that 

Y + 1 

(-1) rIV ^ is also positive definite. 



The system (1) is stable if II is a stability matrix; 
that is if, and only if, all eigenvalues of H have negative 
real parts. 
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The details of the derivation of (5) and (6) are given 

in Appendix A. The solution to (6) can be obtained by 

iteratively solving the r+1 sets of n(n+l)/2 simultaneous 

2 

linear equations or by a technique suggested in Ref. 3. 

The details of this technique are given in Appendix B. 

Evaluation of the performance measure (2) using equa- 
tions (5) and (6) eliminates the need for determining the 
explicit time response of the system but requires knowledge 
of the initial condition, s_(0) . Since a priori knowledge 
of this vector is not generally available, the question 
arises as to how to pick s_(0) in order to arrive at a mean- 
ingful value for J. 

One technique, suggested in Ref. 5, is to determine the 
average value of J over a linearly independent set of ini- 
tial condition vectors. This is equivalent to assuming 
that the initial condition vector is uniformly distributed 
on the surface of a unit hypersphere and results in a de- 
sign which is optimal in an average sense. Another tech- 
nique, suggested in Ref. 6, is to determine the largest 
possible J for every s_(0) within or on a hypersphere of 
radius p and select the configuration which yields the 
smallest of those maximum values. This results in a design 
which is optimal in a worst case sense. For the latter 
scheme, it is shown in Appendix C that for a given design 
parameter set {P}, the largest value for J is given by 
where 



2 

Since the matrices 
only n(n+l)/2 elements of 1 



, V , , . . . ,V are each 
each need to be solved 



symmetric , 
for. 
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/ 



J 1 ({P » VWC-l) r+1 r!W< P »] (?) 

and where ^ max [‘] represents the largest eigenvalue of [•]. 
Here, the functional dependence of J 1 upon the design 
parameter set {P} is stressed. 

D. MINIMIZING THE QUADRATIC PERFORMANCE MEASURE 

To this point we have assumed that the design parameter 
set {P} is unbounded. In general, however, we must allow 
for the possibility of a constrained {P}. There may be 
explicit as well as implicit constraints. The explicit con- 
straints are those which are the result of the physical 
limitations on the system components while the implicit 
constraints are those which result from restrictions placed 
on system performance by the designer 

In either case, the constraining relationships must be 
included in the problem formulation. One means of accom- 
plishing this is to add a large penalty to the cost of a 
given configuration if that configuration requires the vio- 
lation of any of the constraints. With this in mind we de- 
fine the optimal parameter set {P } as that parameter set 
which minimizes J^({P}) and the minimum value of J ({P*}) 
as J ; that is 

J = J, C ( P >) = min J ({P}) . 

{P} 1 

Therefore"^ 



Since r! is a constant in any given design specification 
it docs not contribute to the optimization of the parameter 
set and can therefore be deleted from the performance measure 
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J = min{A 

' {P> 



max 



[(-i) r * 1 v r , 1 ({p))]). 



( 8 ) 



This indicates that the following iterative procedure 
can be used to find {P } once r and Q have been specified. 

1. Select an initial estimate for {P} insuring that 
any constraints imposed are satisfied and that H is a sta- 
bility matrix. 

2. Using equation (6) or the method described in Ap- 
pendix B, find 



select a new set of parameters {P}. If the constraints are 
not satisfied, add a penalty to the cost and select a new 
parameter set. 

5. Repeat steps two through four until some criterion 
for ending the search is satisfied. 



4 

For the examples used in this thesis, the subroutine 
"DIRECT" (NPS Computer Facility) that performs the pattern 
search method of Hooke and Jeeves was used. 



C-i) r+ 1 v r+ i({P». 



3. Evaluate 
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III. OBSERVER DESIGN 



In the design of feedback control systems it is often 
desirable to have information about plant states which is 
not normally available through direct measurement. It is 
sometimes possible to construct a system called an observer 
which can estimate a linear transformation of the plant 
state and, if the transformation is invertible, it is then 
possible to reconstruct the plant state vector itself. 

Traditionally, observer design has been accomplished 
through the use of techniques such as those described in 
Section IIA. In this section, techniques for the design of 
full-order and reduced-order observers by using the quad- 
ratic performance measure as the design specification are 
presented . 



A. FULL-ORDER OBSERVER DESIGN 

We consider the time-invariant, linear system 

x(t) = Ax(t) + Bu (t) 

X(t) = Cx(t) 



( 9 ) 



where x(t) is the nxl state vector, y_(t) is the mxl output 
vector, u(t) is an £xl control vector, A is an nxn constant 
matrix, B is an nx£ constant matrix, and C is an mxn con- 
stant matrix. We consider also the full-order identity 
observer shown in Fig. 1. The observer state equations as 
developed in Appendix D are given by 

£(t) = (A-GC)z(t) + G^(t) + Bu(t) CIO) 
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Figure 1. Full-Order Observer. 



where z_(t) is the nxl observer state vector, £ is the nxm 
observer gain matrix (in this case the elements of G con- 
stitute the design parameter set {P}) , and where it is as- 
sumed that (£,A) is observable. 

Defining the error state as 

A 

e(t) = z(t) - x(t) 



and differentiating with respect to time yields 

e(t) = z^(t) - x(t) . 

Substituting (9) and (10) for x(t) and z_(t) , respectively, 

* 

and solving for _e(t) yields 

e(t) = (A-GC)e (t) . 

We define the matrix F as 

A 

F = (A-GC) (11) 



which yields 



£(t) = Fe_(t) . 



( 12 ) 



We define the quadratic performance measure as 



oo 

A f _r 1 
= 1 t e_ 



(t) £e ( t) dt . 



(13) 



The system (12) and the performance measure (13) have the 
same form as (1) and (2), respectively. Equation (13) can 
therefore be rewritten as 

J = (-l) r+1 rle T (0)V r+1 e(0) 
where, by using (6), 
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T 

F v 



s s 1,2 , . . . ,r 



F^i + VjF = -£. 

Furthermore, the algorithm described in Section IID can be 

* 

applied to find £ , the observer gain matrix which is opti- 
mal with respect to the performance measure (13). 



B. REDUCED -ORDER OBSERVER DESIGN 



We again consider the time- invariant , linear system 
given by 



x(t) = Ax(t) + Bu (t) 
l(t) = Cx(t). 



(14) 



We define the vector x(t) as 

l(t)' 
w(t) 



x(t) = 



where w(t) is an (n-m)xl vector given by 



w(t) = Dx(t) (15) 

and where D is an (n-m)xn matrix, arbitrarily selected so 
that 

-1 

C 

_ D 

exists. By proper selection of the state variables C will 
have the form 



c = [ i : o ]. 

A notationally convenient choice for D is 

d = [ o : i_ ] . 
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We can find w(t) , the estimate for w(t) by constructing the 
reduced-order observer as shown in Fig. 2. For the above 
choice for C and D, the observer state equations developed 
in Appendix D can be written as 

iCt) = (A 22 -GA 12 )z(t) + [A 21 -GA 11 + (A 22 -GA 12 )G] Z (t) 



+ (B 2 -GB 1 )u(t) 

where 

w(t) = _z(t) + G^(t) 



and where the matrices A and B are partitioned as 



A = 



rll. ^ .ri? 

—21 T —22 



B = 




(16) 



(17) 



Here, A^ is an mxm matrix, A^ 2 is an mx(n-m) matrix, A_ 21 
is an (n-m)xm matrix, A 22 is an (n-m)x(n-m) matrix, B^ is 
an mx£ matrix, and B 2 is an (n-m)xJl matrix. Therefore, 
_x(t) , the estimate for _x(t) (which is, for the case 
C = [ I_ ! £ ] and D = [ £ : _I ] , the estimate for x(t)) , 
becomes 



x(t) = 

Defining the 

e(t) = 

and substituting 
tively, yields 



I(t) 

i(t) + G^(t) 

/\ 

error in the estimate w(t) as 
w(t) - w(t) 

(15) and (17) for w(t) and w(t) , 



respec- 



e_( t ) = z_(t) + (GC-D) x (t) . 
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Figure 2. Reduced-Order Observer. 



Differentiating with, respect to time, substituting (14) and 
(16) for x(t) and .z(t) , respectively, and solving for e_(t) 
yields 

e.(t) = CA 22 -GA 12 )z(t) + [A 21 -GA n + (A 22 -GA 12 )G]y(t) 



+ (B 2 -GB 1 )u(t) + (GC-D) (Ax(t) + Bu(t)). 
We note that for C = [ ] an ^ D = [ 0_ : X 1 



GCAx(t) - ( G A ^ ^ • GA^ 2 
-DAx(t) = ["A_ 2 ^ : - A_ 22 
GCBu(t) = GB^uCt) 
- DBu (t) = -B 2 u(t) . 



]x(t) 

]x(t) 



We note also that for y(t) = Cx (t) 



A 21 y(t) = [A 21 : 0]x(t) 
■£A n y(t) = [-GA n : 0]x(t) . 



These relationships, when substituted into the above expres- 
sion for e_(t) , yield 



e(t) = (A 22 -GA 12 )z(t) + (A 22 -GA 12 )GCx(t) 

- [ 0 : (A 22 -GA 12 ) ]x(t) . 



We note that the last term in this expression can be written 
as 



- [ 0 : (A 22 -GA 12 )]x(t) = - (A 22 -GA 12 )Dx(t). 

Substituting this relationship into the expression for £(t) 
and combining terms yields 

e(t) = (A 22 -GA 12 )(z(t) + (GC-D)x(t)) 
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where 



z(t) + (GC-D)x(t) = e(t). 



Therefore e_(t) becomes 



f -.( b ) ^—2 2 — 12 ^ ^ ^ * 



If we define 




( 18 ) 



then 



e(t) = Fe (t) . 



(19) 



Defining the quadratic performance measure as 



CO 




( 20 ) 



o 



we see that the system (19) and the performance measure (20) 
have the same form as (1) and (2), respectively, and that 
(20) can therefore be written as 



Once again, the algorithm described in Section IID can be 
* 

used to find G , the observer gain matrix which is optimal 
with respect to the performance measure (20). 

C. EXAMPLE OBSERVER DESIGN PROBLEM 

To illustrate the application of the algorithm to ob- 
server design we consider the plant whose transfer function 



J = (-l) r+1 r!e T C0)V r+1 e(0). 



By using (6) we have 




s = 1,2 



F T Vi + V-jF = -Q. 
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is given by 



Y(S) _ 1 

0(51 S 3 + 12S 2 + 47S + 60 



The state and output equations for this plant can be written 
as 



— 








— 


0 


1 


0 




0 


0 


0 


1 


x(t) + 


0 


-60 


-47 


-12 




1 










L — 



x(t) = 



£(t) = [100 ]x(t) . 



u(t) 



Here , 



> 










— — 




0 


1 


0 




0 


A = 


0 


0 


1 


» B = 


0 




-60 


-47 


-12 




1 












_ 



, C = [ 1 0 0 ]. 



It can be easily shown that (C,A) is completely observ- 
able. A second-order observer with arbitrary eigenvalues 
can therefore be constructed to provide estimates for the 
unmeasured states C t ) and Xj(t). These estimates, x„(t) 

A 

and Xj(t), are given by 



x 2 (t) 

x 3 (t) 



= w(t ) = z_(t) + Gy_(t) 



We require an initial estimate for the observer gain 
matrix G which yields a stable observer. The observer F 
matrix is given by (18) as 



- A 22 ~ GA.12' 



(18) 



Here 



22 



-22 



0 1 
-47 -12 



a 12 « [ 1 0 ] , 




Therefore 



F = 



‘Si 



■47-g. 
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where the observer eigenvalues are the eigenvalues of £ and 
are given by the solution to 



X 2 + (12+ gl )X + (12 gl +g 2 +47) = 0. 



We arbitrarily select the initial estimate for the ob- 
server eigenvalues to be -7, -2. The characteristic equa- 
tion is therefore 



(X+7) (X+2) . 0 



or 

X 2 + 9X + 14 = 0. 

Equating coefficients of the like powers of X and solving 
for the £ matrix yields 




A parameter in observer design which must be limited is the 
observer speed since, as the observer gains are adjusted to 
cause the errors to die out more rapidly, the observer be- 
comes more sensitive to noise. One measure of the observer 
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speed is the sum of its eigenvalues. It can be easily shown 
that this sum is numerically equal to the trace (the sum of 
the diagonal elements) of the F matrix. 

As our constraint on observer speed we therefore set a 
limit, T, to the most negative value that we will allow the 
trace of the matrix F to assume. To observe the effect of 
T on the resulting design we will allow T to take on the 
values -9, -12, and -15. We see that for our initial esti- 
mate for £ that the trace of £ is -9 and that this £ will 
therefore satisfy any of the values of T previously stated. 

The design algorithm of Section IID was used to compute 

* 

the optimal observer gain matrix G , the corresponding value 

* 

of the performance measure J (£ ) , the observer eigenvalues 

and an ^ the worst case initial condition error vector 

for various choices of the design specifications £ and r in 

the performance measure (20) and for various choices of the 

constraint T. The plant and observer were then simulated 

on the IBM 360 digital computer for each design specifica- 

* 

tion using the corresponding £ and the worst case initial 
condition error. 

Table I shows the effect of varying the design parameter 
* * 

r on £ , J (£ ) , and the observer eigenvalues A^ and A£. For 
these cases the matrix £ is fixed at 




and T is fixed at -12. Note that in each case, the value 
* 

for g^ is zero. Since the trace £ must be greater than or 
equal to T we see that 
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- 12 - gl > - 12 



is satisfied only for less than or equal to zero. The 
* 

fact th.at is always zero indicates that the constraint 
is always active. 

The effect of varying the parameter r on the observer 
error time response is shown in Fig. 3. Note that as the 
value of r is increased, the estimate error for both states 
reaches zero sooner. For the estimate error for the state 
Xj(t), however, this is at the expense of larger error 
values in the early part of the operating interval. 

In Table II the effect of the design parameter on 

* it 

G , J (G ) , X^ and X^ is shown. Here, the value of r is 
fixed at zero and the value of T is fixed at -12. The ma- 
trix has the form 



a - 



q n 0 

0 1 



Again we note that the constraint T is active. 

The observer error time response is shown in Fig. 4. 

The effect of increasing q^ on the plant states is similar 

to that of increasing r, although somewhat less dramatic. 

In Table III the effect of the design parameter on 

* * 

G , J (G ), X^, and X^ is shown. Again, the value of r is 
fixed at zero and the value of T is fixed at -12. The ma- 
trix Q has the form 



^ - 



1 0 
0 q 22 
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Table I. The effect of r on G , J (G ), X^ and X£. 

2 - [o i] • T - - 12 



*11 


1 


2 1 


10 


* 

Si 


0.0 


0.0 


0.0 


* 

82 


-35.0 


-30.0 


-8.9 


J(G*) 


1.04 


1. 50 


3.58 


A 1 


-10.9 


-10.4 


-6.0-jl. 


X 2 


- 1.1 


- 1.6 


-6.0+jl. 


Table II. 1 


rhe effect 


of on 


G*, J(G*), 




r = C 


1, T = -12 
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Figure 3. Observer Estimate Error Time Response as a Function 
of Time Multiplier Exponent r. 
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Figure 4. Observer Estimate Error Time Response as a Function 

of q ir 
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Once again we see that the constraint T is active. 

The observer error time response is shown in Fig. 5. 

The effect of q 2 £ on the response is quite different from 
the effects of r and q^ as we might expect. By placing 
greater emphasis on the cost associated with states which 
represent higher derivatives we do not allow them to devi- 
ate as far from zero and thus tend to slow the convergence 
process down. Hence, as q 2 £ increases, the observer error 
states both reach zero at later values of time. 

It should be pointed out that all observer configura- 
tions were simulated for a period of two seconds for the 
sake of uniformity in the graphical presentations. If the 
response of Fig. 5 were continued for a greater length of 
time, the curves would all eventually approach zero. 

In Table IV the effect of the design constraint T on 
* * 

G , J (G ) , X i and X^ is shown. Again, r is fixed at zero 
and Q is fixed at 

"l o' 

a - 

_0 1 

* 

By noting the values for g^ we see that the constraint T is 
always active. 

The observer error time response is shown in Fig. 6. 

We see that the effect of making T more negative is to im- 
prove the shape of the response at the earlier values of 
time. Note that in the case where T has the value -9 that 
the observer error in state ^(t) actually gets worse for 
early values of time. 
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Table III. The effect of q ^2 on £ > J (G ), X^ and X^. 
r = 0, T = -12 



* 

*1 

* 

*2 
J (G*) 



-3.0 

27.1 

8.30 

4 . 5- j 4 . 2 
• 4 . 5 + j 4/2 



-12 

0.0 

-35.0 
1.04 
-10.9 
- 1.1 



-15 

3.0 

-44.9 
0.16 
- 11.8 
- 3.2 



Table IV, 



The effect of T on G , J (G ) , X^ and X^, 
1 O' 



Or 



o 1 



r = 0. 
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Figure 5. Observer Estimate Error Time Response as a Function 
of q 22’ 
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KEY 





Figure 6. Observer Estimate Error Time Response as a Function 
of Constraint T. 
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IV. OBSERVER/ CONTROLLER DESIGN 



A. TRADITIONAL DESIGN TECHNIQUE 

The traditional method for the design of feedback con- 
trol systems utilizing observer estimates of unmeasured 
states is to specify desired plant performance in the manner 
described in Section IIA as if all states are available 
through direct measurements. A feedback network or controller 
is then constructed which causes the closed-loop system to 
satisfy these performance specifications. The controller 
has available as inputs either the observer estimate of the 
entire state vector or the measured plant states and the ob- 
server estimate of the remaining states. The observer is 
designed to be stable and to have a characteristic equation 
with eigenvalues which satisfy the somewhat nebulous require- 
ment that the observer by "slightly faster than the plant" 

[ 2 ]. 

In this section, an alternative method is described in 
which the observer and controller are designed simultaneously 
using a performance measure similar in form to (2) in Sec- 
tion IIB. This results in a configuration in which both the 
controller design and the observer design are optimal with 
respect to the specified performance measure. 

B. FULL -ORDER OBSERVER/CONTROLLER DESIGN 

The full-order observer for the linear, time- invariant 
system 
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( 21 ) 



x(t) = Ax(t) + Bu(t) 

L(t) « Cx(t) 
is, as shown in Appendix D, 

z(t) = (A-GC)z(t) + G^(t) + Bu(t) . (22) 

We define the control law as 

u(t) = Kx (t) 

where K is a constant £xn matrix. 

It can be shown [2] that the stability of a closed-loop 
system using an observer output as a substitute for the 
plant states in the control law is unaffected by a stable 
observer. The control can therefore be approximated by the 
control estimate, u(t) , where 

u(t) = Kz_(t) . (23) 

The structure for this full-order observer/controller is 
shown in Fig. 7. 

Defining the error state, e^(t) , as 

e(t) ^ z(t) - x(t) (24) 

and differentiating with respect to time yields 

e_ ( t ) = 2 ^(t) - x (t) . 

Substituting (21) and (22) for x(t) and £(t) , respectively, 
and solving for £(t) yields 

e(t) = (A-GC)e(t). (25) 

Substituting (23) into (21) yields 

x ( t) = Ax (t) + BKz(t) . (26) 
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X 




Figure 7. Full-Order Observer/Controller. 
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Solving (24) for z_(t) and substituting into (26) yields 

x(t) = (A+BK)x(t) + BKe(t). (27) 



Defining the augmented state vector, £(t) , as 



5(t) = 



x(t) 
e (t) 



and using (25) and (27), it is easily shown that 



i(t) 



A+BK BK 

0 A-GC 



i(t) • 



(28) 



The performance measure can be defined as 



~ (29) 

J = J t r [x T (t)Qx(t) + e T (t)Re(t)]dt. 
o 

This performance measure is chosen because it accommodates 
the specification of convergence properties of both the 
closed-loop system states and the observer states and im- 
plicitly limits the expenditure of control effort. Equa- 
tion (29) can be written in terms of £(t) as 



J 



J tf(«) 

o 



a a 

0 R 



£(t) dt . 



(30) 



Equations (28) and (30) are in the same form as (1) and 
(2), respectively. It has been shown [2] that the matrix 
T, defined as 

A+BK BK 
0 A-GC 
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has as its eigenvalues, the eigenvalues of A+BK and the 
eigenvalues of A-G£. The matrix Y is therefore a stability 
matrix if both A+BK and A-GC are stability matrices. It 
can be shown that the matrix 

a a 

0 R 

is positive definite if both Q. and R are positive definite 
[7]. Therefore (30) can be rewritten as 

J = C-l) r+1 r!£ T (0)V r+1 i(0) (31) 

\ 

where, by using (6) 



A+BK 


BK “ 


T 


"A+BK 


BK “ 


0 


A-GC 


V + V 

-s + 1 -s + 1 


0 


A-GC 


L J 


• 


c 


- 



A+bk 


BK ' 


T 




'a+bk 


BK 




"a 


o' 


0 


A-GC 


Li 


* 'Ll 


0 


A-GC 




0 


R 



Furthermore, the algorithm described in Section IID can 

* * 

be applied to find K and £ , the feedback and observer 
gain matrices which are optimal with respect to the per- 
formance measure (29). 

A slight modification to the design algorithm causes a 
significant reduction in the number of computations re- 
quired on each iteration. To see why, we note that if the 
observer estimate initial condition vector is assumed to be 
the null vector, that is 
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1(0) = 


0 


then, from (23) , e_(0) 


becomes 


e(0) = 


-x(0) 



and the augmented state vector £(0) becomes 



1(0) = 


x (0) 

(32) 

■x (0) 



Partitioning the symmetric matrix V ^ as 



— r+l * 


hi —12 

T 

V V 

-12 -2 2 


where each element is 


an nxn matrix, then substituting (32) 



and (33) into (31) yields 

(34) 



J = 


(-l) r+1 r!x T (0) [V 11 -V 12 -V 12 T +V 22 ]x(0) . 


We now define the nxn 


matrix V' - as 
—r+l 


- 


[V u -V 12 -V 12 T + V 22 ]. (35) 



Therefore 



J = 


(-l) r+1 r:x T (0)V; +1 x(0). (36) 



The matrix has n eigenvalues, while the matrix Y. t +i 

has 2n eigenvalues. We see, therefore, that the effort 
required in finding can be significantly reduced by using 
(35) and (36) in step three of the design procedure. 
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C. REDUCED-ORDER OBSERVER/CONTROLLER DESIGN 
For the linear, time -invariant system 



x(t) = Ax(t) + Bu(t) 
Z.Ct) = Cx(t) , 

we define the control law 

u(t) = Kx(t) 

where K is a constant Hxn matrix. 



( 37 ) 



(38) 

We define the vector 



x(t) 



A 



L(t) 

w(t) 



where w(t) is an (n-m)xl vector. We let 

w(t) = Doc (t) 

where J) is an (n-m)xn matrix, selected so that 

_ -.-i 

C 

D 



(39) 



exists. By proper selection of the state variables, £ will 
have the form 



C = [1:0]. 

A notationally convenient choice for D is 

D = [0:1 ] . 



A 

For this choice of C and D, w(t) , the estimate for w(t) , 
can be obtained by constructing a reduced-order observer, 
as developed in Appendix D, whose state equations are given 
by 



39 



Z(t) - {^22 ^- 12 ^—^^ + 2 1 "^^1 1 + (— 2 2 ” — — 1 2 ^ ^ X.(b) 



+ (B 2 -GBi)uCt) 



( 40 ) 



where 

w(t) = z_(t) + G^(t) . (41) 

The matrix K and the vector x(t) can be partitioned as 

K = [K m :K°] , x(t) = 

where, for the stated choice for £ and D, 



x m (t) 

x°(t) 



£(t) = w(t) = x°(t), w(t) = x°(t). 



The control law, given by (38), can therefore be approximated 

b X _ 



u(t) 



K 



x m (t) 

x°(t) 



or 



u(t) = [K m :K°] 

or 

u(t) = K m y(t) + K°z_(t) + K°Gy(t) . (42) 

The structure for this reduced-order observer/controller is 
shown in Fig. 8. 

We define the observation error, £(t) , as 

e(t) ^ w(t) - w(t) . (43) 

Substituting (41) for w(t) and (39) for w(t) yields 



___Z_(t) 

z(t) + Gx(t) 



e_(t) = _z(t) + Gy(t) - Dx(t). (44) 
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Figure 8. Reduced-Order Observer/Controller. 



I 



I 

» 

f 
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Substituting (37) for £(t) and rearranging yields 

e(t) = z(t) + (GC-D)x(t). 

Differentiating with respect to time yields 

e (t) = £(t) + (GC-D)x(t) . 

Substituting (37) and (40) for x(t) and ^(t) , respectively, 

» 

and solving for £(t) yields, as shown in Section IIIB, 

e(t) = (A 22 -GA 12 )e(t) . 

Solving (44) for z^(t) + G^(t) and substituting into (42) 
yields 

u(t) = K m Z (t) + K°(e(t) + Dx (t) ) . 

We note, however, that 

Dx(t) = x°(t) 

and 

X.(t) = x m (t) . 

Therefore , 

u(t) = Kx (t) + K°e(t) . (46) 

Substituting this estimate for u(t) into (37) yields 

x(t) = Ax(t) + B(Kx(t) + K°e(t) 
which, after rearranging, becomes 

x(t) = (A+BK)x(t) + BK°c(t)). (47) 

We define the augmented state vector, C(t) , as 



Using this definition and the expressions for x(t) and _e(t) 
given by (47) and (45), respectively, yields 

rO 



, £(t) - 



A+BK 

0 



Br 



- 22"— 12 



£(t) . 



(48) 



The performance measure can be defined as 



(49) 



J = J ~ t r [x^ (t) Qx (t) + e^ (t) Re (t) ] dt 



o 



which can be rewritten in terms of the augmented state vec- 
tor, j[(t) , as 

'"a o' 



= / t r d 



(t) 



0 R 



I(t)dt. 



(50) 



Equations (48) and (50) are in the same form as (1) and 
(2), respectively. Furthermore, the matrix 

A+BK BK° 

- - 22"— 12 

is a stability matrix if both A+BK and ^ 2’— 12 are 
matrices. Also, the matrix 



a o 

0 R 



is positive definite if both Q and R are positive definite. 
The performance measure (50) can therefore be rewritten 



where, by using (6), 



A+BK 


BK° 


T 








‘a+bk 


BK° 




0 


-22"— 12. 




Vi 


+ 


V s*l 


0_ 


-22"— 12_ 












• 

• 

• 








s = 


"A+BK 


BK° 


T 








"a+bk 


BK° 




0 


—2 2" —1 2 _ 




Si 


+ 




0 


— 22 "— 12_ 


■ 



V 
— s 



1 , 2 ,. 

"a a 

0 R 



Furthermore, the design procedure of Section IID can be 

* ife 

applied to determine K and , the feedback and observer 
gain matrices which are optimal with respect to the speci- 
fied performance measure (49) . 

Once again, it is possible to reduce the number of 
computations required for each iteration of the design 
algorithm. In this case, we set the observer estimate, 
w(0) = 0_, which yields £(0) = -Dx(0) . Therefore ji(0) be- 



comes 



ICO) = 



x(0) 

-Dx (0) 



Substituting this into (51) yields 



J 



(-D r+1 r! 



'x(O) 

-Dx (0) 



V 



r+1 



(52) 

'x(O) " 

-Dx(0) 



Partitioning the (2n-m)x(2n-m) symmetric matrix V ^ as 



Vi 



t— i 
i—i 

>i 


—12 


v -u 


y .22 
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where is an nxn matrix, V^ 2 an nx ( n ' m ) matrix, and 

V _22 is an (n-m)x(n-m) matrix and substituting into (52) 
yields, after expanding 



(53) 

J = (-l) r+1 r!x T (0)[V n -D\ 1 /-V 12 D+D T V 22 D]x(0). 



We define the nxn matrix V' as 

— r+1 

ZN-1 = [V n -D T V 12 T - V 12 D ♦ D T V 22 D], (54) 

Therefore, (53) can be rewritten as 

J= (-l) r+1 r!x T (0)V; +1 x(0). (55) 

The matrix has n eigenvalues, while the matrix Y. r +i 

has 2n-m eigenvalues. We see, therefore that the effort 
required in determining in step three of the design pro- 
cedure can be reduced by using (54) and (55). 

D. EXAMPLE OBSERVER/CONTROLLER DESIGN PROBLEM 

To illustrate the application of the design algorithm 
to observer/controller design we consider the plant whose 
transfer function is given by 



Y ( S) _ 1 

0T5T S 2 + 15S+50 ‘ 



The state and output equations for this plant can be written 
as 



x(t) 





0 1 ~ 




V 


1! 




^c(t) + 






-50 -15 




1 



u(t) 



£(t) = [10 ]x(t) . 

Here , 
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, C = [ 1 0 ]. 



Then 



A = 



A+BK = 





' 0 


1 “ 






Y 


= 






, B = 






-50 


-15 






i 












[k x 




" 0 


1 ’ 




0 




= 






+ 








-50 


-15 




_ 1 _ 












_ 





0 1 



|k 1 ~ 50 k 2 -15 



Also 



where 



F A 22 £Ai2 



A 22 = [-15], A 12 = [1] , G = [g] 



Therefore 



F = - [ 15 + g] . 



We must select an initial estimate for both K and £ 
which causes the matrices A+BK and £ to be stability matrices. 
As in the case of the observer alone we must constrain sys- 
tem speed. In this case, however, we sum the eigenvalues 
of the closed loop system. It is easily shown that this 
sum is given by the trace of the matrix 

"A+BK BK° 



—22 "—1 2 



which is given by 

k 2 - 30 - g. 
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We arbitrarily select the trace constraint T to have a value 
of -25. Therefore 



-25 < - 30 - g. 

We see that the open-loop system is stable. We therefore 
choose an initial estimate for K of [0 0] . We must now 
choose a value for g which satisfies both 



and 



-25 < - 30 - g 
-(15+g) < 0 . 



Or, combining the two inequalities, 



-15 < g < - 5. 

We choose for our estimate of G, therefore, 

G = [-10] . 

The design algorithm of Section IID was used to compute 

* 

the optimal controller feedback gains K and observer gains 

it 

G , the corresponding value for the performance measure 

it it 

J (K ), the plant and observer eigenvalues A^, and A_, 
and the worst case plant and observer estimate error ini- 
tial conditions for various choices of Q and r in the per- 
formance measure (49). The constraint T was fixed at -25 
and the design parameter R was fixed at 

R = [1] 

in all cases. The plant and observer were simulated on the 
IBM 360 digital computer for each performance measure spe- 

it it 

cified using K , G , and the corresponding plant and esti- 
mate worst case initial conditions. 
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Table V shows the effect of varying the design parameter 
* * * * 

r on K , G , J(K ,G ), A^, A and A^. In this case the £ 
matrix is also fixed at 




0 

1 



Figure 9 shows the effect of varying the design parameter r 

on the closed-loop system response to the worst case initial 

conditions. We see that increasing r improves the system 

settling time in both states at the expense of larger values 

of t ^ ie earl y part of the operating interval. 

Table VI shows the effect of varying the design parameter 
* * * * 

q^ on K , G , J (K ,£ ), A^, A^, and A^. In this case r is 
fixed at zero and £ has the form 



2 = 




0" 

l 



Figure 10 shows the effect of varying q.^ on the closed-loop 
system response to the worst case initial conditions. We 
see that- as the value of q^ is increased, the settling time 
improves for both states at the expense of larger values of 
x^(t) in the early part of the operating interval. 

Table VII shows the effect of varying the design parame- 

^ A A * 

ter q £2 on K , £ , J (K ,£ ), A^, A^, and A^. In this case 
r is fixed at zero and £ has the form 



2 = 



1 0 
0 q 22 
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r 


0 


1 


2 


* 

k i 


25.7 


-17.1 


-217.4 


* 

k 2 


- 9.4 


- 9.0 


- 8.0 


* 

g 


-14.4 


-14.0 


-13.0 


* 

,G ) 


1.02 


0.34 


0.06 


h 


- 0.6 


- 1.0 


- 2.0 


X 2 


-23.3 


-20.8 


-11.5-jll 


X 3 


- 1.0 


- 3.2 


-11.5+jll 



Table V. The effect of r on K , G , J(K ,G ), X^, 

'l o' 



a = 



0 1 



, T = -25, R = [1] 



and X 



3’ 



q ll 


1 


2 


10 


* 

^1 


I 25.7 


-15.5 


-25.9 


* 

k 2 


- 9.4 


- 9.5 


- 9.0 


* 

g 


-14.4 


-14.5 


-14.0 


* 

> £ ] 


1.02 


1.46 


3.37 


x i 


- 0.6 


-0.5 


- 1.0 


X 2 


-23.3 


-23.0 


-20.3 


X 3 


- 1.0 


- 1.5 


- 3.7 



Table VI. The effect of q.. on K , G , J(K ,G ), X 1 , X~, 
and X 3> r = 0, T = -25, R = [1]. 1 Z 
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Figure 9. Closed -Loop System Time Response as a Function of 
Time Multiplier Exponent r. 
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x 2 





(b) Closed-Loop System State x^Ct). 



Figure 10. 



Closed-Loop System Time Response as a Function 
°£ q n . 
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Table VII. The effect of q 97 on K ,G ,J(K ,G ), X n , X 9 , 
and X 3 . r=0 , zz T=-25, R=[lJ 1 z 
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Figure 11 shows 
system response 
see that as q 2 £ 
gish, requiring 



the effect of varying q^ on the closed-loop 
to the worst case initial conditions. We 
is increased, the system becomes more slug- 
longer settling times, as we would expect. 
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I 



KEY 




Figure 11. Closed-Loop System Time Response as a Function 
of q 22 . 
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V. CONCLUSIONS 



A new technique for observer and observer/controller 
design which requires the minimization of an integral of 
time-multiplied-square-error (ITSE) performance measure has 
been developed. An algorithm which is well suited to digi- 
tal computation and which is especially useful for higher- 
order systems has been presented. An example of a computer 
program which implements the algorithm for the case of an 
observer/controller is shown in Appendix E. The computa- 
tions are performed off-line with a typical computation time 
of thirty seconds. 

It has been demonstrated through the use of examples 
that modification of the time response of an observer or 
a closed-loop system utilizing an observer/controller to 
initial condition error can be achieved by a suitable change 
in the appropriate design parameters of the ITSE performance 
measure. Furthermore, it has been demonstrated that the 
qualitative change in the system time response reflects 
directly the qualitative change in the performance measure 
design parameter. 

It is recognized that for second-order and, perhaps, 
third-order single - input , single-output systems the rela- 
tionships between traditional design specifications and 
system transfer functions arc simple and relatively easy to 
use. In these cases, the design procedure suggested above 
is generally not efficient. However, when dealing with 
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higher-order, multiple-input, multiple -output systems, the ad- 
vantages of the new technique become more significant. For 
these systems, the design specification and system transfer 
function relationships are highly complex and the traditional 
methods are therefore extremely difficult to use. However, 
the new technique offers a relatively straightforward ap- 
proach to design which does not become more complex as the 
system order, the number of inputs, or the number of outputs 
increases . 
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APPENDIX A 



00 

EVALUATION OF THE INFINITE TIME INTEGRAL f t r £(t)Qs (t) dt 

o 



For the system of linear, first-order differential 
equations 

s(t) = Hs(t) (1A) 

if there exists a scalar function 

VjCsCt)) = sj (t)V 1 ^(t) (2A) 

subject to the conditions 

1) V 1 (0) = 0 

2) V^(s_(t)) > 0, for all s_(t) f 0_ 

3) V 1 (s(t))^ «>, as ||s(t)|| -*■ °° 

4) V^(s_(t)) < 0, for all s_(t) f 0^ 

5) V^sCt)) = 0, s (t) » 0 

then the system (1A) is asymptotically stable [4] . 

Since V^(s_(t)) is a scalar, it is equal to its transpose. 
Hence, from (2A) , we see that the matrix is symmetric. 

We now consider 

VjCsCt)) = A^tjv^sCt) + ^(tjv^Ct). 

Substituting (1A) for £(t) yields 

^(sCt)) = ^(tjH^sCt) + s^tjV-jHsCt) 

which, after rearranging terms becomes 
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V x Cs Ct) ) = s T (t) [h t v x + V^lsCt). 



We now define the matrix Q as 



q = - [H 1 ^ + VjH] 



and note that since is symmetric, then Q must also be 
symmetric. Therefore 



At this point we note that equation (2A) insures that 
conditions one and three are satisfied. Condition two is 
satisfied if is positive definite and conditions four 
and five are satisfied if Q is positive definite. 

Integrating the left hand side of (3A) over the finite 
time interval, T, yields 



VjCsCt)) = -s T (t)£s(t) . 



(3A) 



T 



T 




o 



o 



T 




o 



If (1A) is asymptotically stable, then 



T->oo 



lim V 1 (s(T)) = 0. 



Therefore , 



oo 




O 



5 8‘ 



Or, from (3A) , 



00 

j s T (t)2s (t)dt = s_ T ( 03 V 1 s_ 



(0) 



(4A) 



where is the solution to the algebraic matrix equation 



HV 1 + V^H = . 



(5A) 



We now consider 



gf[t£(t)V 1 £(t)] = i T (t)v 1 £(t) + t [£ T (t) V^Ct) + s T (t)V 1 s_(t) ] 



= v 1 (s_Ct)) + tv^^CsCt)) 



Hence 



CO OO 00 

J 3 ^[ts T Ct)V 1 s (t) ] = J V 1 (s(t))dt + J tV 1 (s(t))dt 

o oo 

It can be shown that if s_(t) is the state vector of an 
asymptotically stable, free system, then for a finite n, 



OO 

J d[t n s T (t)V 1 s_(t) ] = t n s T (t)V 1 s_(t) 



= 0 



Therefor e 



OO 

■hi 



(5(t))dt + 



OO 

/ tv l<i 



(s(tD)dt 



or 



j" tV 1 (s(t))dt = -J V 1 (s(t))dt 
o o 
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or, substituting (3A) for V^(s_(t)), 



oo oo 

f ts T (t)2s(tjdt = - J VjCsCOJdt. 
o o 



Now suppose there exists a scalar V 2 (s_(t)) such that 

V 2 (s(t)) £ s T (t)V 2 s(t) 



and 



V 2 (s(t)) = V 1 (s(t)). 



Then 



or 



V 2 (s_(t)) = s_ T (t)V 2 _s(t) + £ T (t)V 2 s 
V 2 (s(t)) = s T (t)[H T V 2 + V 2 H]s(t). 



However, from (8A) and (2A) we have 

V 2 (s(t)) = s T (t)V 1 _s(t) . 

Therefore, we see that 

!! T v 2 * v 2 h = 

where, since is symmetric, then V ? is symmetric. 
(8A) we have 



CO CO 

- J V 1 (s(t))dt = - J V 2 (s(t))dt 

o o 

co co 

- / vpsCOJdt = - J dv 2 (s(t)) 

0 O 



(6A) 



(7A) 



(8A) 



From 
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oo 

- f V 1 (s(t))dt = V 2 (s(0)) 
o 

oo 

- J" V 1 (s(t))dt = s T (0)V 2 s(0) 



Substituting (6A) yields 



OO 

■J ts T (t)C>s_(t)dt = s T (0)V 2 s(0) (9A) 



where V 2 is the solution to the algebraic matrix equations 



H V- 


+ V 0 H 


= V, 


— —2 


— 2— 


—1 


T 


+ V^H 


= - £• 



(10A) 



The argument can be continued for higher powers of t. 
Thus for a stable system (1A) and a symmetric, positive 
definite Q, 



OO 

f t r s T (t)gs(t)dt = (-l) r + 1 r!s T (0)V r+1 s(0) (11A) 



o 



where is the solution to the simultaneous algebraic 

matrix equations 



H V A1 + V , 1 H = V 
s + 1 — s + 1— — s 



H Vl + W H = - q. 



s = 1 , 2 , . . . , r 
(12A) 
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APPENDIX B 



A METHOD FOR THE DIRECT SOLUTION OF THE SYMMETRIC MATRIX V ^ 

— r+1 

A procedure for generating from the system of 

simultaneous linear algebraic matrix equations 

^♦l - V^H = ^ s = l,2,...,r 

h t v x + V^H = - 

is described in Reference 3. The procedure, which is es- 
sentially repeated here with minor amplification and nota- 
tional changes involves the direct solution for V ^ and 
bypasses the need for solving the (r+1) (n+l)n/2 simultaneous 
equations . 

We first form the matrix S_ as follows: 

1) Form the symmetric array of integers 

1 2 4 7 

2 3 5 8 

T = 4 5 6 9 

” ss 7 8 9 10 

’n(n+l)/2 

where n is the dimension of H, Q, and V , 

2) Form the n intermediate square matrices, S ^ , 
i=l,2,...,n whose elements arc equal to one and whose order 
is n(n+l)/2. 
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3) Proceed along the first row of I_ ss and for each 
integer, j, that is missing from the sequence 

{1,2 , . . . ,n(n+l)/2} , set all entries in the row and column 

of equal to zero. 

4) Proceed row by row through the H and matrices 
simultaneously. For the first non-zero entry in the first 
row of , enter the value h^, for the second non-zero 
entry in the first row of , enter h^* and so forth until 
all entries of H have been used to replace all the non-zero 
entries of S^. 

5) Multiply the first column of by two. 

6) Proceed along the second row of I. as in step three, 

making the appropriate entries in the matrix. 

7) Repeat step four for the S_£ matrix. 

8) Repeat step five for the second column of the 

matrix. 

9) Continue the procedure until all of the matrices 
have been generated. 

10) Form the S matrix from 



S = 




We next form the vector v , , from 

— r+1 



!r.l * i' (r+1) a 



where v r+ ^ and are the leading diagonal elements of the 
symmetric v r+ ^ and Q matrices, respectively. That is, 



6 3 



v ll" 




q n 


V 1 2 




q 12 


V 2 2 




q 2 2 


V 13 


a = 


q 13 


V 2 3 




q 2 3 


V 33 




q 33 


>_ 




q nn 



The matrix V r+ ^ can now be generated from the vector v_ r+ ^ 
in a simple manner. 
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APPENDIX C 



EVALUATION OF THE QUADRATIC FORM s T (0) (-l) r+1 r!V r+1 s (0) 

Let x be an n dimensional, unit length vector. That is 

x^x = 1. (1C) 

The n dimensional vector y_ whose direction is the same as 
that of x but whose length, p, is arbitrary is therefore 
given by 

= px . C 2 C ) 

We consider the quadratic form 

C = £ T A£ (3C) 

where A is symmetric. Substituting (2C) for y_ yields 

C = p 2 x T Ax . (4C) 

Let us constrain to lie on or inside an n dimensional 
hypersphere of radius R. That is, 

0 < p < R. (SC) 

We assume we have no further a priori knowledge of the 
vector y_ and therefore wish to consider its worst case value, 
that is the value of y_ which causes C to take on its maximum 
value for a given A. 

We sec from (4C) that C is maximized for a given A when 

T 

p is a maximum and when the scalar quantity x Ax assumes its 
maximum value. From (5C) we sec that the maximum value for 
p is R. Therefore, 
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(6C) 



C = R 2 max{x T Ax}. 
max — — 

We define the quantity, q as 

q = max{x^Ax} (7C) 

It can be shown [8] that the matrix A in the form 

x^Ax (8C) 

T T 

where x x = 1 and A = A, has n real eigenvalues, A^ and n 
corresponding real eigenvectors x ^ , i = l,2,...,n. Further- 
more, the maximum value of the form (8C) is the largest 
eigenvalue of A. That is 

q = max{x T Ax} = * max [A]. (9C) 

Also, the vector ’which gives this maximum value is the eigen- 
vector corresponding to the largest eigenvalue. 

Substituting (9C) into (6C) yields 

C max - RZ \»*W- O0C) 

Thus, C is the value for C in (3C) which represents 
the worst case value for y_ for a given A. We wish to mini- 
mize C as a function of A for this worst case condition. 

Since R is a fixed quantity, it plays no role in the mini- 
mization of C and we can therefore assume that it is equal 
to one. This is equivalent to assuming that y_ lies on the 
unit hypersphere. 

We are concerned with minimizing the quantity 

J = s T (0)(-l) r+1 r!V r+1 s(0) (11C) 
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r +l 

where the matrix (-1) r!V r+ ^ is symmetric and positive 
definite and where s^(O) represents the worst case initial 
condition of the state vector s_(t) . We can assume, without 
loss of generality, that s^(O) lies on the unit hypersphere. 
From the argument above, we see that J^, the maximum value 
for J for a given parameter set {P}, is given by 

J = X [ (-l) r+1 r ! V Al ] (12C) 

1 max 1 '- J — r+l J v J 

and that the worst case initial condition vector, s^(0) , is 

r +l 

the eigenvector which corresponds to r '— r+1^ * 
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APPENDIX D 



OBSERVER THEORY 
A. THE NEED FOR OBSERVERS 

Given a free plant characterized by the state and output 
equations 

x(t) * Ax(t) (ID) 

X.(t) = Cx(t) (2D) 

where x(t) is the nxl state vector, y;(t) is the mxl measure- 
ment vector, A is a known nxn matrix, and C is a known mxn 
matrix, we wish to know the value of the state x(t) • If C 
were an nxn matrix representing n linearly independent mea- 
surements, then C ^ would exist and x(t) could be found from 

x(t) = C'^t). (3D) 

This is not usually the case, however, so we will assume 
that there are fewer linearly independent measurements than 
there are states. That is, we will assume that m is less 
than n. 

Since C ^ does not exist for m less than n, an alterna- 
tive scheme must be found to determine x(t) . One method of 
estimating the state vector is to construct a model of the 
plant and operate it under the same conditions as the plant. 
Such a model is called an observer. The following sections 
develop the state equations for two different types of 
observers . 
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B. FULL -ORDER OBSERVERS 



Let the observer state equations be given by 



z (t) = Az (t) . 



(4D) 



The solution to (ID) is given by 

x(t) = e— t x(0) 



(5D) 



and the solution to (4D) is given by 



z (t) = e— z (0) . 



(6D) 



Therefore, if z_C 0 ) = x(0) , the observer state will be iden- 
tical to the plant state for all t greater than or equal to 
zero. Since, however, in general ^(0) does not equal x(0) , 
the identical dynamics of (ID) and (4D) would result in the 
persistence of the initial condition error and the observer 
state would not provide an accurate estimate of the plant 
state . 

An alternative observer is one which contains not only 
the plant model but also a feedback forcing function made 
up from an error signal. Since the only accurate indica- 
tion we have of the states is the measurement vector y;(t) , 
let us assume that the error signal represents the difference 
between the plant measurement vector and an identical 
"measurement" of the observer state vector. Specifically 
we assume that the model and the feedback error signal are 
represented by the observer state equations 



The feedback is provided by the nxm observer gain matrix G. 
Rewriting (7D) in terms of x(t) and z(t) we see that 



z(t) = Az ( t) + G( Z (t) - Cz(t)) 



(7D) 
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z(t) = Az(t) + GC(x(t) - _z(t)). (8D) 

Defining the error state e^(t) as 

e(t) ^ z(t) - x(t) (9D) 

and substituting into (8D) yields 

z(t) = Az(t) - GCe(t) . (10D) 



It is clear, therefore, that by using the proper choice for 
G, the error, e_C t ) , can be used to drive z_(t) toward x(t) . 

More precisely, differentiating (9D) with respect to 
time yields 

e(t) = z.(t) - x(t) . 



Substituting (ID) and (10D) for x(t) and z_(t) , respectively, 
yields 

£(t) = A£(t) - GCe (t) - Ax (t) 

£(t) = A(_z(t) - x(t)) - GCe (t) 
e (t) = (A-GC)e (t) . (11D) 



That is 

e(t) = Fe(t) 
where we define as 

F = A - GC. 

The solution to (12D) is given by 



(12D) 



(13D) 



e(t) = e-e(O). (14D) 

We see, therefore, that if the eigenvalues of £ all have 
negative real parts, then for a finite e (0) , 
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lim e_( t) = 0. 

t->oo 

It is shown in Ref. 4 that the eigenvalues of the ob- 
server can be positioned arbitrarily if the pair (C,A) is 
completely observable. That is if the matrix 



nr * t t • 

[ c 1 : a 1 c 1 : 



(A T j‘C 



T: • • • : (A T ) n_ 1 c T ] 



has rank n. From (13D) we see that the eigenvalues of F are 
controlled by the selection of the observer gain matrix (j. 

The above argument is easily extended to a plant which is 
driven by the input u(t) . The plant state and output equa- 
tions become 

x(t) = Ax(t) + Bu(t) (15D) 

£(t) = Cx(t) (16D) 



and the observer state equations become / 

z(t) = Az(t) + G(y(t) - Cz(t)) + Bu(t) (17D) 

where u_(t) is an £xl control vector and | is a known nx£ 
matrix. 



c. reduced-ordf:r observers 

The full-order observer possesses a certain amount of 
redundancy since we estimate some states which are essentially 
being measured directly. We therefore desire to construct 
an observer which estimates only those (n-m) states which 
are not being measured directly. We start with the state 
and output equations 

x(t) = Ax(t) + Bu(t) (19D) 
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£(t) = Cx(t) . 



(20D) 



We define the new state vector x(t) so that its first m 
components are the measurement vector £(t) and the remain- 
ing (n-m) components are the vector w(t) . That is 



x(t) 



4 



Z(t) 

w(t) 



(21D) 



If w(t) is selected so it is some linear transformation of 
the state x(t) , namely 



w(t) = Dx(t) . 



(22D) 



Then, in view of (20D) and (22D) , (21D) can be written as 



x(t) = 




x(t) . 



Furthermore, if D is chosen such that 



(2 3D) 



exists, where 




C ' 

D 

then, the state x(t) can be found from x(t) . Namely 

x(t) = M ^x(t) . 

Assuming that M ^ exists, the task becomes one of obtain- 
ing an estimate for w(t) . Writing the state equations for 
x(t) we have 
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£(t) = Ckc(t) = CAx ( t ) + CBu(t) 
w(t) = Dx(t) = DAx (t) + DBu(t) # 

That is 

x(t) = MAx(t) + MBu(t) . 
Substituting (25D) for x(t) yields 

|(t) = MAM” 1 x ( t) + MBu(t). 



(26D) 



(27D) 



(28D) 



Defining 



A' = MAM -1 



B' = MB 



then x = A’x + B'u(t) 

which can be written in partitioned matrix form as 



(29D) 



let) 


II 


w(t) 









A • 
A ll 

A ' 
A 21 



A ' 

. 1 ? 

A ' 
A 2 2 



LCt) 

w(t) 



M. 

—2 



u(t). 

(30D) 



Therefore 

£(t) = A’^Ct) + A’ 2 w(t) + B’u(t) C 3 ID) 

and 

w(t) = A^y^t) + A£ 2 w(t) + B^u(t). (32D) 

Since y_(t) is a measured quantity, its time derivative, 
y_(t) , can be determined. The input u(t) can also be con- 
sidered as a known quantity. Therefore, solving (31D) for 
A-j^wCt) yields 



Ai 2 w(t) = £(t) - A' lZ (t) - BJu(t) 
where the elements on the right represent "known" quantities. 
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The quantity A^ 2 w(t) can. therefore be used as the measurement 
vector for a full-order observer for the plant whose state 
equations are given by (32D). We therefore write the ob- 
server state equations as 

w(t) = A^wft) + AJ^Ct) + BJu(t) 

( 33D) 

+ GCy(t)-Aj LlL (t)-Bj L uCt)-A« 2 ^Ct)) 



where (3 is an (n-m)xm matrix, w(t) is the estimate of w(t) , 
and where the first three terms represent the plant model 
and the last term represents the feedback term. Rewriting 
(33D) yields 



w(t) 



(A« 2 -GA* 2 )w(t) + (A^-GA^lCt) 
+ (B’-GBpu(t) + G£(t). 



C 34D) 



It has been shown [4] that if C C_, A) is observable then 
(A^ 2 ,A£ 2 ) is also observable. Therefore G, the observer 
gain matrix, can be selected so that (A 22 -GA^ 2 ) has arbitrary 
eigenvalues . 

Since differentiation is usually troublesome to perform, 
it is desirable to eliminate the term y_(t) from (34D) . To 
accomplish this we subtract the term Gy(t) from both sides 
of (34D) to get 

w(t)-Gi(t) = (A’ 2 -GA’ 2 )w(t) + (A^-GA'^xCt) 

+ (B£-GB_i)u(t) . ( 35D) 



Adding and subtracting 
of (35D) yields 



(A 22 -GA^ 2 )G}l(t) from the right side 
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w(t)-G£(t) = (A] 2 -GA^)(w(t)-G Z (t)) 

+ (± 21-^11 + CAJ 2 -GAi 2 )G )}1 (t) 

+ (B' -GBpu(t) . (36D) 

We now define 

z (t) = w(t) - G^(t) . 

Equation (36D) can therefore be rewritten as 

z(t) = (A' 2 -GA* 2 )z(t) + (A^ 1 -GA' 1 +(A' 2 -GA^ 2 )G) Z (t) 

+ (B’-GBJ)u(t) (37D) 

where w(t) can be computed from 

w(t) = z(t) + G^(t). (38D) 
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APPENDIX E 



SAMPLE COMPUTER PROGRAM FOR OBSERVER/CONTROLLER DESIGN 

A computer program is given on the following pages which 
can be used to design an observer/controller in the manner 
prescribed by the algorithm of Section IID. It is assumed 
that the plant state and output equations are given by 

x(t) = Ax(t) + Bu(t) 

£(t) = Cx(t) (IE) 

and that the observer state and estimate equations are given 

by 

z(t) = (A-GC)^(t) + G£(t) + Bu(t). (2E) 

The control law is given by 

u(t) = Kz(t) (3E) 

and the performance measure to be minimized is given by 

oo 

J = jf t r [x^ (t)Qx(t) + e^ (t) Re (t) ] dt . (4E) 
o 

The following symbols are explained in the order of their 
appearance in the input section of the program: 

BETA - The initial step length for the minimizing search 
routine DDIRCT. 

EPS - The minimum step length allowed. 

D - The minimum value that the sum of the eigenvalues 
(trace) of the closed- loop system is allowed to take on. 
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ALPH - The weighting multiplier which is used to assign 
a cost to any violation of the trace constraint. 

N - The order of the matrix A in (IE). 

JY - The maximum number of searches the minimizing search 
routine is allowed to perform. 

NR - The time multiplier exponent r in the performance 
measure (4E) . 

M - The row dimension of £ in (IE) . For full-order 
observer/controller design M is set to zero. 

A - The matrix A in (IE) . 

Q - The matrix £ in (4E) . 

BB - The matrix B in (IE) . 

GG - A vector whose elements are given by the initial 
estimate of G and K in the form 

' G " 

K 

SS - The R matrix in (4E) . 
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HIS program implements the design algorithm of section iid 

OR THE OBSERVER/CONTROLLER. 
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SUBROUTINE PRINTM IA,l,M) 
DIMENSION A ( 7 1 7 ) 

DOUBLE PRECISION A 



oooooooooooo 

«-< CM m la O CD O' O ^ AJ 






•» 




«— * 


p-H 




««» 


II 




LA 


“5 




• 


•» 




r\j 


« — » 




«— 4 






UJ 


•» 




•» 






X 






ro 


< 










CO 






•» 


— * 




«i 


o<\J 






•» •» 






^4vOLU 




w 


II — z> 






— ^ 


2 


K- 


’JJ^ 


a 


< 


— <1— h- 


z> 


5: 


— 2 : 


i— 


QdO 


o cr o 


LU 




020 


QC 


LL lU 






r\j 



o o o 



101 



oooooooooooo 

^CsirO^iAsOh-oo^O*— *c\j 



-J 








9% 








< 






in 








• 


< 






CM 


> 






*— 4 


K- ^ 


— -» 




LU 


e o 






•* 




<W> 




X 


0£r*-00 


< 




m 








•* 


<o 






* 


LU LU 


— J(\J 






zza: 


•» •» 




- 


hQQ. 


< vO 'JJ 




— * 


I— 


II —3 






3oO UJ 




2 


i— 


0^3 


UJ ■— < 


CrC 


< 


CtUJCO 


—< h - h- 


3 


r: 


COSI3 


— e 


>— 


rx o 


3^U 


oq:q 


UJ 




00 oo 


020 


cc 


LL UJ 




H 




CNJ 



o o o 



102 



DOUBLE PRECISION VERSION OF THE SUBROUTINE OIRECT 



000000000000000000000000000000000000000000000003 
in 3 <0^0 3 (\icn st in vO 300^0— in 3 od O' o—«c\irn^ in 'O 3 coo>o— nicov* in so 300 
— r-< — -h — — — — »-i c\i <\j c\j cnj cm c\j c\i r\j cm cnj [T\ rn o rn rn m m m m rn n* st sfr st sfr <t* st ^ n± 









2 














3 










a — 


in 




1 1 






03 


— 


3 






1 


UJ 


CNJ 




— Ql 






O 03 




2 






1— 1 


2 UJ 


— 




23 






1- 03 


3 


OO 






21— *0 


— O 


• 


2 


— 00 






<< 


3 


2-12 


to 




n<Qh 


< O 




< 


X 




UJ 


3 > 


Ql 


i— <32 LL 


LU 


2 


LU 


3 <1 


LU 


X 


LU 


X 


o 


O 


3 


0X3 0 


O 


2 


a *2 - 


1- > 


3 


3 


LU CD 


O 


a. 


• < 03 


3 


233 


CO 


•» 


UJ2— 00 


< 


O 




CO 


cc 


X 


XX — > 


O 


O 33 


< 


> 


JO< 


O 


< 


oO 


1— 


LL 


< 


— 3 3 — 




3X<3 


i— i 


UJ 


Jm| — •» 




> 


00 


ooo 




X 


• XOO 1— 


O 


OOOCL O 


2 


X 


— I — 1—2 


— 00 o 




LU 


H-3 


h- 


LU 


>— oCO 


►— i 


3 < 1 < 


< 


< 


LL3<0 


OO LU — 


UJ 


O 


X 


►— < 




023 33 


3 


233> 


> 




-J — 


302 


X 




•* 


X 


cO 


3 — 0 00 3 


00 


2 00 




*> 


00 O 00 1 — 


wOQLU 


3 


00 


— • LU 


UJ 


O 


l-X —33 


O 


0200 


2 00 


o 


— 00 3 O 


oo<X 




— 


MX 




LL 


< 00 • 3 


2 


03 3 C 


UJ 


a: 


02 


— 3 


• 




x< 


2 




— 1—2300 


3 


3 2 


UL > 


> 


3U IO 


322 




LU 


Ql2 


O 


•> 


000322 


< 


030— 


OLU 


2 


mQ< LL 


00 C 


— 




w 


Ql 


Q 


32 — 30 C 




2 < 2 3 


UJ 


a 


> 


LU > LU 


• 


— 


00 LU 


3 


JJ 


X h< — 


a 


< — 




2 


• LU 3 


-J 00 


3 


00 


X 




1— 


X 00 <03 •* 




-OX 


00 


*> 


0033O 


003 


o 






O' 


O 


— — 3— < O 


2 


3231- 


o 


00 


LUC OO 


<o 


• 


CL 


• 2 


UJ 


3 


033 * 


►— t 


3 — 0 — 


*2 


*> 


— J X LU 00 


X<3 


o 


UJ 03 LL 


00 


fr- 


3— <203 


< 


0002 


2< 


o 


CO — C03 


oo a 


X 


1— LU— O 


3 


ill 


Q 00 > — < O 


3 


< — 2 


o 


wJ 


<1- -J 


22 


a: 


00 33 




a 


<0.3 > • 


00 


>QC H ► 


►— ( UJ 


3 


*— <oooo aQ 


oo a 


• 




COX 


>- 




XOO 00 3 > 


o 


O w fr— 


32 


■JJ 


232< 


ZLLQ 


3 


1 — 22 <1 


CO 


00 


•» w3— 3 




2 


uO 


O 


< — .HH 


— 3 


o 


2 — 3 2 




< 


—00 O 2 X 


o 


O 3 2 X 




•• 


>3<2 


<LU 


• 


UJZILO 


a 




• — • oo 0< 


3 


— X2 3 


ox 


o 


oo3< 


3 2 • 


o 


22 O 


LU 


o 


3 2C^00»— <2 




3 33 


3 w 


X 


i-22;» 


2<3 


w 


2o— 12 


1— 


Cu • 


023333 


Q 


O 1-3 




2 


2 — 0 


Ooou 




31— <10- 


00 


03 


•QCQC0OO3 


3 


2 33 < 


<X 


»• 


33 03 


O— LU 


2 


O 


23 


UJ 


QCO 


033X02— 


00 


3<< 


o 


Ql 


O 2 


2 


O 




oo cC 3 


1- 


a: • 


X33333 


• 3 


3 CIO 


32 


< 


2X1- LU 


3X — 


I— 02— LUOO 




3X 


aC 3 0 2 00 3 *Q 


3 2 


o< 


o 


33- O 


01-0 


O -U LU 


1— • 


LUO 


QL X — X33 


3 O X < 


UJ 


-j 


a-hh 2 


LUO 


<hIXX2^Q^ft: 


— 3 || 3 3 00 3 


o<o 


Xoo 


LU 


3333 


22> 


LL00 2 


OLU OlLJ 


1 LULUaC 


-a: *0X3 2 


x< 


3 


o 


O 03 


— OcO 




iU 


c£ H00G>l-OOG .00(3— — 


021- UJ — 


23 


*• 


002a: LULU 


QO 


23 • 


<21-31—23 


• O<aco0a:203 


3oo 2 


— o 


— 


2 — 320 


O 


002 UJC O 


ox 


3 *2>Ql3 — ! 


CD 33 — 20 


2 UJ 


00 


LU 00 •— 1 2 a 2 2 LU 


— 3000 2322< 


003 2oO cQ X 


2 — 


3 3 — 


►— < 2 


3 


h^DQMiuaiiJct: 


1— 00 — 


U.DIO 


arm 


• • -J> O X 


OO 


3 00 3 3 


2 — 


00 


LU i 


njq:i— LU 


o 


23031— 1— 3< — 


Ql3 *3223033 *0 3 3 < • 


o 


•» 


X3>X3- llooi— 


300LU x 


C — Ql •COOX — 22i 


— < >* 


O < 2 O 


o 


2 


COcOOOX O 


0 — 1 — h'JJ'tVL) 


i—o *o< ►< - 20 a: 


•'CO 0 2 


C LL 


*> 


a: a: 


— 1 — oc 


LU 


►— * 


X CQ — • 


C000320I-XQ 


— < — < *-o> 


OO 


— t 


< 2 >- 3 2 2 — < 


2in2 


»•< • 


o>- 


• 


3X3 *3320000 


| —332 


O 


00 


Cl 3 0 LU— X— 00 


CM0020020 


-L-< 2 — 0200 


• 2>XO 


OO 


CL 


1— Jl- CD 233 — 


• Ql 


• 


O tUUJ 


1 — UJ r-M 


OO — 3 m — 3 3 < O 3 2 


o 


w 


LL o <r — • 2 


— 1— uu 


O J'JJNh 


i ct 


1 


II • X 


oo>x 


<x 


3 


03— x O LUI— 2 CL 001— a: 2 LUX — J 2 — 


II 


<2323 


• 


3 


o 


>1— LUX 


ao< — < — oooujqi— 23 <0 




X — i 


0<2 


LU UJ 


LU 


2 


OO 
















1— 2 


2 


o 1 1 


1 1 o 


| 


1 


1 


1 




1 


1 




< 


►— ( 


— 


30 
















OLU 


Q 


3 


Ql Q2 






O 










LU O X 




X 


< * 


— < 


o 




a: 




> 






l/l Jl- 


O 


>— < 


— o • • 


O 


o 




> 




3 






o 


O 


2 — 


00 0 3 03 


o 




2 




X 






o-o> 


LUC 


OOO 


OlLU *<X 


LU 




o 




< 


2 




| — CO 


OO 


oo a. 2 


loq 2 02 


O 


00 


2 




X 


2 




o 


< 


LU 




















CL 


00 

3 


a 





















oooooooooooooooooooooooooooooooooooooooooooooooo 



103 



oooooooo OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
O cm m <4* in vO r- co o> o <m m -t tA o co cn o —* c\j m sj* la no h» oo o o — < cm m 4- in o r- <x> O' o -h cm ro vf • n o 
>*• la la la la in la in la la in o o vo o 'O no sO so 'O h- r- h- h- h- I s * r- r- oo oo oo oo oo oo <x> oo oo oo O' o o o> cr> o\ o> 



X O t 


< 
















LU OO'O 






2 












O -J 


2 




X 












X <r 


O 




•» 












LU X) 






> 














O 


00 


LU 












r» Q 


LU 


H 


X 












»— 1 LU 


h- 


►— t 


< 












i 002 


< 




2 












t X 


2 


X 


•* 












UJ UJ2 


O 


LU 


O 












-J 1-0 


>— « 


O 


X 












• • LU 


00 


2 


> 












h-2 2oO 


LU 


3 


2 












ox < — 


O 




Q 












UJ (X 




2 


X 












xxi < 11 — 


• LU 


< 


•» 












— 00 ft-O 


LA CO 


X 


00 












Q< CL 


H 


to 


»- 








•» 




o LL 1 — 


J- 


oo 


O 








rH 




>• OO 


00 oo 


LUX 


-J 








II 




COUJ o 


-» O 


XCL 


-J 








l-H 




X h- 


2 


N 


LU 








•» 




Oh- 


00 


3oo 


O 








M 




LU O Q — J 


LU h- 


QfX 


•» 








— 




X2LUI-< 


-J M 


LULU 


o 








•» 




00 — h- 3 


co 


Xoo 


X 








2 




-« h-O • 


< •* 


3 


X 












JOOM oo 


»-* a 


00 


•* 








*• 




x<2 *-2 


X or 


22 


X 








-J 




s: ooh 


< < 


<I-H 


< 








> 




O LU «0 


> a 


X 


o — 








LU 




LJOOJLUUJ 


2 


Oo 


-J CO 








X 




o<X2co 


h- < 


o< 


LU CM 












<«*.< • 


2 h- 


X2 


O — 








2X 




21 


LU «N Si 


XX 


•* o 








*• « 




00 oo 00 X o 


OO 


COLU 


HH U 








Oh 




mmLU X 


2 LUh— 


3b- 


00 00 








-J II 




>oo< 


UJOOO 


00 X 


X 








O—i 




I— J— OUJlu 


0-02 


UJ 


00 


o 






UJ *• 




oos: oo oo 


UJCL 


2 




CvJ 






O — 




0.0. < 


05 00 


oo 


X*— oo 












M-LUOUJ 


2 — — 


»— MJJ 


•*rM(M 


O 






o— 




33X ct: 


i— i 


l-X 


*-H | — 


h- 






IH 




00020 


LUsO 


o< 


OOQh 


o> 






xoo 




-J— «LL 


u-x 


2-1 


0- -X 


CvJ O 






►X 




O **X UJ 


0<2 


30 


— xx 


CM O 






X w 




2 — tX • cQ 


3 • 


LL LU 


1 


O - nJ* 






< •* 




• LU • 


X oo — O 


O 


b- <T •* 


h— Q> •“ — * •* vQ 






o^ 




**LU LULU 


LU h- O X 


o 


o— — 


cvirOOO •* 






-JOO 




OO h-O 


CO 2UJ< 


2LU • 


X co X 


O •* ^QfM LA 


X 




LUX 




II • o< 


5h2:o 


<C0LU 


*— » ft — 


O— •*> •* 


< 




QOO 




25 


3< 


2 


O-J — 


f^CMCM • ^ UJ 1 ^ 


O 




00 






2X1— — 1 


oo< 


o <r lo__ i 


*-*CT> (C\J X O 




>— i 








1-30 


LUh- 2 


UJX< 


IAAI— .O' • C — O 


UJ 


00 


cma>o 






200 0.01 


2 


LUX > 


r-H *»X CM LU «—■ . 2 _Jin 


XO 


X 


•*roro 






3Zhh 


— X 


2 2 LU 


• o< oo > 






00 ^ 


—4 




002 03 2 


h— *~h LU 


M 1— O 


hfMO- • -J II UJ || 


— < II 


oO ^ 


vQvQ 


00 




^Huoa 


3'^q. 


V— — — X 


O J 30 J X 


II 






X 




XX o 


0^0 


3 O OO LU 


• ^tUIIUj<f J 




II II 


2 


00 




«OlL 


X— X 


0 — 20 


X^QXXO>2> 


M 




X UJ UJ 






552m 


COOOX 


X — 1 UJ LU 


1 1 





— -J 


*— 1— h- 


II 




LU 


3 


tQQ.5h 


X X 


o 


00< 








X 


00 


32 — 2 


LL UL U. LL LL LL < U_< 


Q-J 


x> 


XXX 


00 








00 — O — 


———^^—2—2 


OOO 


00 UJ 




00 










— ICMrn LA 


sOh- 




00 


a> 



XoO 

•*X 



II II 



Oh 

Hw 
►— » 
OX 
O X 



ooooooooooooooooo o 



oo o o 



o o o 



104 



oooooooo oooooooo oooo ooo oooooooo oooooooooooooooo 
k oo o> o cm tn so k oo o o —i c\i ro >4" m so I s - oo c> o ph cm o 4- in so h- oo c> o —1 c\i co -4- m o h- oo o' o ph c\i .t* 
O' ^ ^ O O O O o O O O O O ^ ^ ^ ^ ^ ^ ^ ^ ^ cnj c\i cnj cnj cnj cvj r\j rn rn rn rn ro rr, rn rn rn cn vt nJ- 

H HH rH *— < iH H H H pH f— < rH f— I *— < pH f— 1 1 — ( *-"4 < H rH *— < pH pH H pH H pH pH pH I pH pH pH rH pH H pH pHHpH pH »— H pH i—4 



o o 

_l -J 

00 00 



I I 









II 


II 














•» 


o 


•H 


p^* 














pH 


ro 


HH 


HH 














II 




-*p 


H-P K 














HH 


O 


o 


O LU 














•* 


P"p K 


-J 


-J X 




o 


pH 








-p* 


CM 


00 


oo H 




CM 


CM 








HH 


O 




1 














— p 


OO 


— 


— — k 




o 


o 








HH 




sO — 


-p* HH 




K 


K 








X 


p- * 


pH pH 


rH W 














a 


O-J 


. 


M 




M Oh 


o 




— 




W 


o> 


O' *-H 


— X 




o OO 


o 




O 




ro » 


UJ 


CM O 


oo a. 




-1 -1 






CM 


CD 


pH 00 


— X 


**d 


Cl — 


CO 


00 -^00 






H 


M 


h00 


hh<t 


LA • 


• hO 


rH 


+ oo 1 


00 pH 




f-i 




lA 


oox 


pH h* 


h-PH— O 


^ rH 


—►"■Hi h— — * pH 00 K — * 


pH 00 O 


K 


pH 


o 


rH ^ 


Cl t 


o 


jHt-nO 


HH 4-Q 


MMiLLJm + • UU pH 4- • _J 


'XI 


— 


»— 




i/VJJ 


t 


— r • 


X -IK 


X x — 1 K I X — 1 LLJ 00 


X 






c\l ro 


•O 




p-* hh Cl CM 


a. < 


*‘HOKCL< JhCL<0 1 


K 


•* 




pH *» 


UU • 


i— < ^ H 


HH 00 


00^p>-h><A 


HIJ -p> • hp 


> • PH 






pH O' 


o 


oo 


ii -pp— r*- 


HP CL II II 


ooooxuucm 


II CL 00 II 00UJ-H || 00UUpH || XCM || UU 


CQ 


pH 


- — . — 


• <r 


^-h O^-hpH 


HH 


CL 


pH X 


X 0l CM 


X 


HH 




X 


oo ;> 


-IX 


X II PHP-s 


II II 00 II 2 


II || || || CL— || 


II CL — O0 


— X 




Oo 


uu 


oo uu 


CO 00 0.00. HH 


o 


CM hh 00hh 


00 PH Qmm 


o 


pH >— 


— 1— 






Spf-HSP 


MM || _J HH 


CMhH-Hj— H 


-J^p— II K 


NPh 


K 


OO 


HH 






UU*hm 


00 X <00 


UUU)HI< phI< O 






00 o 


u-a: 


LL LL 


O'XULOLLXOOX 


CLCLOO>00 


OXJICL>lLXCL>lLU^Oia 


o 


<o 




PH pH 


Ohhhh(J 


HH \ Q_ Q_ 


oooooouu < 


O K 00 CL 00 UU ph O*. 00 UJ hh 00 00 O CL O 


o 




pH (M 


rO sj- 


in 


ok oo 




O' 


O 


pH CM 






pH pH 


pH pH 


pH 


pH pH pH 




pH 


CM 


CM CM 





o o oo 



oo oo 



oo o 



105 



23 IF (KN) 24,25,24 1440 



oooooooooooooooooooocooo oooo oooooooo 
insor^oo(T' 0 ^ovjro<Mn'Or^co^o— 'Ovjrn<MH' 0 ^coo^o— jovjmvtmvor^coo^o 
st >t ^ 't in tn tn in mirun in in in so nO nO no vO no vO vO nO no r*^ I s - h- h* go 

i-H <H i-H w-M «-H ^ i-^ <— ( *H -H i-^ ^ «H r-4H «-H H H«-H ^ «— ( ^ ««*«^ 



o 

h— 

O 

o 



00 



00 

CO 

<t 

o 

■ft 

o 

Q 



X 




in 








•• 




• 


m 






H 




t 








II 




h- 


o 






*— « 




o 


L— 






*• 




• 










O' 




o 






♦— « 






o 






w 


O 


*—* 








►— « 


b- 










X 




»-H 


OCl 






CL 


o 


00 


_i<r 


i—* 




w 


o 


CL 


-JO 






•» 




1 


LU _J 


o 




00 


— . 




QLU 


-1 




'S> 


»— « 


■-H 


• Q 


00 






00 


w 


-ft 


■ft 






o. 


t-M 


JO 


o 




CO 


oo 


^X 


• X 


^x 




m 


• 


•.CL 


CLc£ 


-CL 




•• 


UJ 




< 


■ "4 




vO 


o 


II 00 LU 


O II 


II II 




■— 


• 


-hcQ3 




>-H 


O' 




l/) 


<2 


IJJCL 






UJ 


00 


voa— ■ 


Q< 


00*-» 


O 


1— 




CM— h- 


— O 


CM-^ 


1— 






2 


-J 


O 




CL 


LL 


OlO 


u_aj 


OO 


o 


3L 


I-H 


QmO 


*— «Q 


QUO 


o 




in 


s0 


h- 


CO 




CM 


C\J 


cm 


CM 


CM 







*-2 ► 


— * 




XDQ^I 






NlL<h 


vO 




- >X 


• 




vOX vO 


in 




• in Q ► 






in.— i c£ cm 


LU 






n0 




UJ O •* 


•* 




*■ ►x 


X 




II UJCMm 


s0 




>H- 


CM 




OO t%s0 


w 




ISX ► 


X 




cl mx 


s0 




IO -s0 


• 




mixcM 


in 












— .^.^.UJ 




CM\ CLX 


■ -s0 




»*m <c M* 


sO O 'O 




vOw> •* 


• • *x 




• ► cL 


inminm 




in n Q< 


«• 




^ 2> 


UJLULUr^ 




ULIZI 


sOOO • 




-^vOX 


•* •« *.in 




II -b- 


XXX— 1 


* 


XCMX 


sOsOOUJ 


•» 


cLin«-4v0 


CM CM CM •» 


H 


<j •» #. •» 




II 


OXXCM 


\S» 


HH 


JfMfO*-* 


OsOO 


•• 


LU •• •* •• 


• • •!— 




Q 00 XX 


ininio— < 




x—^m 


hhhX 


<w> 


oo •* 


LULU LUUJ 


»— « 


•« II QCsO 


'O'O'O 


l/) 


X < - 


«*x 


CL 


c\JJ>X 


X X X co 




•*> s0 


mmm^ 


•« 


X LU h- CM 


— — —o 


*— « 


UX^S 




l/) 


CL <£ X-* 


• • • t— « 


CL 


<Xs0X 


intnm •. 


00 


LUX •'<t* 


i— ^ •— ( r»"t 


«. 


00 CO CM *• 


LU UJ LU II 


O 


«<H(^ 




CL 


KX -< 


uj 2 O 


> 


L)(MX> 


2 cL CL CL 


m2 


uj *.m 


— »OUJ> 


mo 


ocso *-x 


001-2 




M *Xh 


*—< CL H— O 


_Jvt 


Qinmx 


odx < 


<rm— 


f—4 1—4 •»xQ 


O'JJCLO 


— » > *•(> 


X’-ULLJ •* 


O O X 


1 o Lumm 


>t *.0£M 


xxxo 


m •» 


— * II _J— t 


CO CO 00 —c 


II CM II C\J II v0 






m m — ^ 


o>x — 



O O 02 2 1 --J I'OXHhhh 

>►->*->*— i-o *- 22 X 2 : 

2 2 2 ►-»►— CC Q— • -cLCL cLCLcCCh 

OOOUOIL^LLI QIhX<OUOOZ 
^O^O^'—'JCCx: LL^-U^>LLU.Li-LL^U 

^CMrOvt 

a> o — <CMm vt m sOh-coo' 

r\J CO ('OrOfOrO m mm mm 



o o 



oo o o o 



o 



106 



ooooooooooo 







GOn^CO 








•'CO •* 








< < 








•w •» 








— 2! 








odt^a 




— * 




h- 




l-H 




szr^oo 




•0 








~ ) 




<0 








UJ UJ 


2 






ZZQC 


•* 


•% 




mQOl 


fH 


— < M 




H-— < 


11 


II 




2)oOUJ 


►— < 




2 


OZJ 






a: 


QC UJcQ 


H 


H *» 


2) 


CD 51 ZD 




►—i 


h-Q 




O 


O — 


UJ^ 


OOOQ 


O 


O CD 


<xlu 



107 



ooooooooooo 
cvj m in vO ^ oo o> o ^ 



N- 



2 — 

♦•I s - 



o 








♦•h— •* 




—5 




CQwcQ 




•* 




•*C0 *• 








< < 








■'-«» •* 




CO 




~2 




4- 




QNO 








<X f*N-l 




—> 




Slf^oo 




•* 




2 -^ 




M 




<tu 




W 




UJ UJ 


2 


2< 




ZZ£^ 




•» 




mOh 




^ II 




V— • 


Ik 


II 




Ol/)UJ 


»— < 




2 


02-1 




O 


oC 


0£ lUCQ 


»-H 


H •• 


O 


CO 2:0 




»— 1 


h-O 


D^O 


O 


O — 


UJ2 


tOOQ 


Q 


OO 


OCLU 



o o 



o 



SUBROUTINE MULT ( A , B , C , L , M , N) 

DIMENSION A ( 7 » 7 ) » B( 7 , 7 ), C( 7 , 7 ), D( 7 , 7 ) 
DOUBLE PRECISION A,B,C,SUM»D 



oooooooooooooooooooooo 
<\j in o r*- oo O' o — ■ cm m >4- in <0 oo ^ o ^ cm 

^ (MCMCM 






CO 







”3 






— •» 


















—1 






Vo 






— 


2 E 








0 


< 








-J 


2:0 


21 + 


CO 


_J 




r* 




•o^ 






► 


rH 


^ • 


■ — 1 ZD 


II 




— < II 


II 


it 0 


II cO 




II 


II 


»— « 




—> 




►— « 






11 


II 








CNJ 


CM 




Vo 


m 


m *■* 




51 




M 1 




►— » 


a 


□ D 


O^) 




0 


O — 


0 


OOO 


OcO 


0 


0 


OO 



«-h f\J m 



00 o o 00 o o 



LU^ 

gclu 



109 



TURN 



LIST OF REFERENCES 



1. Luenberger, D. G., "Observing the State of a Linear Sys- 

tem," IEEE Trans, on Military Electronics , Vol. MIL-8, 
p. 74-80 , April 1964. 

2. Luenberger, D. G., "An Introduction to Observers," IEEE 

Trans, on Automatic Control , Vol. AC-16, p. 596-602, 
December 1971. 

3. MacFarlane, A. G. J. , "The Calculation of Functionals 

of the Time and Frequency Response of a Linear Con- 
stant Coefficient Dynamical System," Quarterly Journal 
of Mechanics and Applied Mathematics, Vol. XVI, p. 
259-271, 1965. 

4. Kalman, R. E., Bertram, J. E., "Control System Analysis 

and Design Via the 'Second Method' of Lyapunov - 1. 
Continuous Time Systems," Journal of Basic Engineering , 
Transactions of the ASME , p. 371-393, June 1960. 

5. Levine, W. S., Athans , M. , "On the Determination of the 

Optimal Constant Output Feedback Gains for Linear 
Multivariable Systems," IEEE Trans, on Automatic Con - 
trol , Vol. AC-15, p. 44-48, February 1970. 

6. Park, C., "Time Domain Design of Observers," MSEE 

Thesis, Naval Postgraduate School, Monterey, California, 
1973. 

7. Kreindler, E., Jameson, A., "Conditions for Nonnegative- 

ness of Partitioned Matrices," IEEE Trans, of Automa - 
tic Control , Vol. AC-17, p. 146"j February 1972. 

8. Bodewig, E., Matrix Calculus, Interscience Publishers, 

Inc. 1959. 



110 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0212 * 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Department Chairman, Code 52 2 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

4. Assoc Professor D. E. Kirk, Code 52Ki 3 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

5. LT Robert Ernest LaRock, USN 1 

SUPSHIP 



Newport News Shipbuilding and Drydock Co. 
Newport News, Virginia 23607 



111 



Thesis 156347 

L269 LaRock 

c.l Optimal design of 

observer/control lers. 



I 4 AUG HO 

250CT 83 

1 ■ OCT 87 



2 5 7 8 7J 

2 8 4 98 

3 3 3 6 9 



Thes i s 

L269 LaRock 

c.l Optimal design of 

observer/control lers . 



15S347 



thesL269 

Optimal design of observer/ controllers. 




3 2768 001 03237 8 

0U0LEY KNOX LIBRARY 



